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Granular materials constitute an intermediate state 
of matter between fluids and solids Jj]. Over the past 
decade, a large series of experimental assays on granular 
systems, concentrated in the quantitative studies of com- 
paction, mixing, segregation and patterns formed by vi- 
bration, have been done [2-9]. They uncovered numerous 
mechanisms (percolation, convection, ordered settling, 
arching, etc.) for segregation of dissimilar grains; re- 
markable standing-wave patterns (stripes, hexagons, dis- 
ordered waves, cluster of localized solitary excitations, 
etc.) in vibrating layers of granular materials and highly 
inhomogeneous and localized stresses chains in nearly 
static granular media. Nevertheless, in contrast to vis- 
cous fingering in fluids, a widely studied pattern-forming 
phenomena [9] , grain-grain displacement seems to be an 
unexplored area in granular physics. Here we investi- 
gate the pattern formation when a grain is displaced by 
another type of grain in a radial Hele-Shaw cell. We 
show that several morphologies can occur, ranging from 
rounded to fingered patterns, interconnected by a con- 
tinuous crossover. Fourier analyses shows that, in con- 
trast to the rounded patterns, the fingered shapes present 
mode selection. 

The experiments were carried out in a radial Hele-Shaw 
cell made of two parallel square glass plates, each with a 
thickness of 0.5 cm and lengths of 90 cm (bottom plate) 
and 80 cm (top plate). The top plate had a hole of 0.5 cm 
in diameter in its centre for the injection of grains. Ini- 
tially, a dense uniform circular monolayer of styrofoam 
spheres was formed on the cell's bottom plate, always 
covering the same area. After that, the cell was closed 
with the top plate, over which weights were placed to 
prevent upward movement, and another kind of grain 
was injected through its hole. In all experiments steel 
spheres with a diameter of 3.94 ± 0.03 mm were man- 



ually injected, whereas three types of styrofoam initial 
monolayers were used: type A, with an average diameter 
of 2.3 ±0.1 mm and low polydispersion; type B, com- 
posed of spheres with diameters in the range of 0.5 to 4.0 
mm and average diameter of 2.9 ± 0.6 mm; and type 
C, with an average diameter of 3.5 ± 0.3 mm and low 
polydispersion. After every 100 g of injected spheres, 
corresponding to the mass of N — 395 spheres, the 
growing patterns were photographed in order to follow 
their evolution. 

Different displacement patterns were observed as a 
function of the diameters of the styrofoam spheres and 
the plates separation. In Fig. 1, a pattern ~ 50 cm wide 
showing well developed fingers, formed after the injec- 
tion of 6715 steel spheres can be seen. In addition to the 
expected interface between the steel/styrofoam spheres, 
there is an additional interface in the displaced grains. 
This second interface is the border between two regions: 
an outer one, still a dense monolayer, and an inner one 
in which the spheres have moved upward and, hence, is 
no longer a monolayer. All the displacement patterns 
generated in the experiments exhibited this second inter- 
face. Also present is a similar upward movement of the 
displacing steel spheres. 

In Fig. 2, a qualitative morphological diagram is pro- 
posed based on the data available from the limited range 
of grains diameters and plates separations D tested in our 
experiments. For large D (> 5.60 mm) the displacement 
patterns are circular. For intermediate D (ranging from 
(5.17 to 5.60 mm), the interface of injected steel spheres 
exhibit small wavy or cusp instabilities for type A (small) 
or type B (intermediate) displaced styrofoam spheres, 
respectively. For small D (< 5.17 mm) fingered pat- 
terns are observed, including fingers having pronounced 
tip splitting, resembling the branched erosion channels 
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formed in water/glycerine saturated Fointanebleau sand 
p0| , for the smallest D. For a given value of D, the fin- 
gers become more apparent as the average diameter of 
the displaced styrofoam spheres is increased. All, these 
different morphological phases are probably not sepa- 
rated by sharp phase transitions and, instead, contin- 
uous crossovers among them might be observed. Indeed, 
in Fig. 3, the patterns formed by the steel spheres in sty- 
rofoam spheres of type B, suggest a continuous crossover 
from a fingered pattern (Fig. 3a) to an almost circular 
one (Fig. 3d), as the plates separation D increases. The 
number of injected spheres (N = 6320) and the scale 
of the drawings are the same for all patterns showed in 
Fig. 3. From now on, the sequence of patterns observed 
in Fig. 3, which represents a cut of constant average di- 
ameter (type B) in the morphological diagram will be 
analysed. Similar behaviours are found for other parallel 
cuts. 

From Fig. 4, in which the distance of each point along 
the border to the mass centre of the pattern is repre- 
sented, one can see the time evolution of the patterns 
depicted in Fig. 3a, d. For the fingered pattern, the for- 
mation and splitting of the fingers can be clearly seen 
(Fig. 4a-c). On the contrary, for the rounded pattern, 
although its average radius increases in time, the devia- 
tion around this value remains small, and fingers are not 
formed (Fig. 4d-f). Figure 5 shows the Fourier trans- 
forms of the displacement borders in Fig. 4. For the fin- 
gered pattern of Fig. 3a, it is possible to see that, as its 
size increases and the fingers become more pronounced, 
modes are selected, which shift from high to low wave 
numbers (Fig. 5a-c). Rauseo et al. (llj have observed 
the same behaviour for the power spectra of the curva- 
ture versus arc length along the interface of the patterns 
generated by fluid-fluid displacement for their regime of 
fast flows. Conversely, for the pattern of Fig. 3d, which 
remains rounded during the growth, no mode selection is 
observed (Fig. 5d-f). 

In fluid-fluid displacement the fingers are formed when 
a low viscosity fluid is injected in a high viscosity one. 
For the systems studied here, we suspect that the role 
of viscosity is played by the friction of the grains against 
themselves and the plates of the cell. The great elastic- 
ity of the styrofoam spheres enhances this friction since 
they can deform and increase their contact area, provid- 
ing a higher "effective viscosity" as compared to the one 
for the steel spheres. This increase in friction becomes 
more critical as the plates separation is decreased, lead- 
ing to more pronounced fingers in the formed patterns. 



In fact, experiments in which the styrofoam spheres were 
replaced by hard plastic beads generated less conspicu- 
ous fingers as compared to the ones formed in styrofoam 
for the same plates separation. The relationship between 
friction and "effective viscosity" can explain why finger 
instability occurs only after a certain number of steel 
spheres have been injected. At the beginning, the resis- 
tance against the displacement front is small since the 
styrofoam spheres, apart from some isolated defects, ini- 
tially form a monolayer. Consequently, its "effective vis- 
cosity" is low (comparable to that for the steel spheres), 
and the border of the displacement pattern is circular 
with small fluctuations superimposed. As the injection 
proceeds, the styrofoam spheres pile up and can jam in 
a crescent number of points, increasing their "effective 
viscosity" . When this value overcomes the "viscosity" of 
the driving grains, the circular interface becomes unsta- 
ble and fingers begin to form. Since the stress generated 
by the injected grains concentrates at the finger tips, the 
jammed styrofoam spheres in these regions are eventually 
released due to their elasticity, leading to a fast growth of 
the fingers. In contrast, if hard plastic beads, which can 
be deformed only slightly, are used, the stress concentra- 
tion at the fingers is not sufficient to loose the jammed 
beads around them, consequently, slowing down the tip 
growth and contributing to stabilise a circular displace- 
ment front. 

The width of the multilayered region depends on the 
plates separation. For large separations the styrofoam 
spheres can easily move upward, releasing the stress 
caused by the displacing grain in a narrow region. On the 
contrary, for small separations the upward movement of 
the styrofoam is hindered and a wider region is necessary 
to release the stress. As the pattern develops, the driv- 
ing front faces an increasing resistance, directly observed 
by the crescent force necessary to inject the steel spheres 
in the cell. At later stages of the experiment, even the 
displacing grains are under a so large stress that the steel 
spheres themselves begin to move upward, initially in a 
region close to the injection point and them spreading 
across the pattern. Differently than for the styrofoam 
spheres, even for the largest plates separation used, the 
steel spheres cannot rise enough to form a bilayer. The 
multilayered region in the styrofoam, as well as in the 
steel spheres, never reaches the outermost ones. There 
is always a monolayer of spheres that seems to be free of 
stress. 

The main result of this experimental study was the dis- 
covery that several morphologies, connected by a contin- 



2 



uous crossover, can be found for grain-grain displacement 
in a radial Helc-Shaw cell. It was possible to observe that 
the formation of fingered patterns proceeds through the 
selection and amplification of characteristic wavelengths. 
We are currently refining the morphological diagram and 
investigating the complex stress distribution in this gran- 
ular system. In particular, the effects of elasticity will be 
the focus of the next experiments. 
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FIG. 1. Fingered pattern formed by the injection of steel spheres in a monolayer of styrofoam spheres of type B. The plates 
separation is 4.93 ± 0.04 mm. 
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FIG. 2. Morphological diagram for the patterns formed by the steel spheres injected in styrofoam particles as a function of 
the diameter of the styrofoam spheres and the plates separation. 
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FIG. 4. The border of the patterns shown in Fig. 3a (a-c) and Fig. 3d (d-f) for three different numbers of injected grains. 
a,d, N = 395. b,e, N = 3160. c,f, N = 5925. 
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FIG. 5. Power spectra of the border of the patterns shown in Fig. 4. a-c, Transforms for Fig. 4a-c. The two strongest 
modes selected correspond to the formation of fingers and their tip splitting, d-f, Transforms for Fig. 4d-f. 
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